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We study a three-dimensional plaquette spin model whose low temperature dynamics is glassy, 
due to localised defects and effective kinetic constraints. While the thermodynamics of this system 
is smooth at all temperatures, we show that coupling it to a second system with a fixed (quenched) 
configuration can lead to a phase transition, at finite coupling. The order parameter is the overlap 
between the copies, and the transition is between phases of low and high overlap. We find critical 
points whose properties are consistent with random-field Ising universality. We analyse the interfa¬ 
cial free energy cost between the high- and low-overlap states that coexist at (and below) the critical 
point, and we use this cost as the basis for a finite-size scaling analysis. We discuss these results in 
the context of mean-field and dynamical facilitation theories of the glass transition. 


Introduction - There are several theories that aim to 
explain the rapid increase in the viscosity of supercooled 
liquids close to their glass transitions m- Some of these 
theories predict that phase transitions should occur, ei¬ 
ther on cooling or in response to some kind of exter¬ 
nal field [US]. Such phase transitions are relevant be¬ 
cause their associated order parameter fluctuations may 
explain the characteristic fluctuations that are observed 
in the dynamics of supercooled liquids [3] . The theory of 
dynamical faciliation (DF) OITU] is based on a dynami¬ 
cal order parameter, and the associated dynamical phase 
transitions occur in ensembles of trajectories iziiiiiiia. 
This theory is encapsulated by a class of simple model 
systems that includes both kinetically constrained mod¬ 
els [T3| and plaquette spin models [14]. A different 
theoretical approach is based on mean-field calculations 
within a replica formalism lilli]. It proposes a static 
order parameter - the overlap - which measures the sim¬ 
ilarity between two randomly chosen configurations of the 
system in equilibrium. This approach predicts a phase 
transition, as found in certain spin-glass models, that oc¬ 
curs when configurations of the system are biased to be 
similar to fixed reference configurations. 

In this paper, we consider a three-dimensional plaque¬ 
tte spin model ~ the square pyramid model (SPyM) [16] 
- whose relaxation behavior is well-described by DF the¬ 
ory. Using computer simulations, we show that this 
model also exhibits phase transitions when biased by 
its overlap to (quenched) reference configurations. The 
properties of the thermodynamic phase transitions that 
we find are consistent with the universality class of the 
random-field Ising model (RFIM), as predicted by field 
theories for generic glassy systems [HI HI]. Evidence for 
this scenario has also been found in atomistic simula¬ 
tions [min]. Thus, at low temperatures, the SPyM ex¬ 
hibits both anomalous dynamic fluctuations (as predicted 
by DF) and anomalous overlap fluctuations (as predicted 
by mean-field). These results are important for two rea¬ 
sons: first, they provide strong evidence that theoretical 
predictions of RFIM criticality [TT] [TS] are applicable in 
three-dimensional systems; second, the fact that a sin¬ 
gle model displays features of both mean-field theories 


and DF theory provides an opportunity for a connection 
between these apparently contradictory scenarios [3]. 

Model - The SPyM [IB] consists of Ising spins Sj = ±I 
on the vertices of a body-centred cubic lattice, with 
periodic boundaries. We use C = {si,...,sn) to in¬ 
dicate a configuration of the system, whose energy is 
^0 = -y where the sum runs over 

upward-pointing square pyramids within the lattice, and 
the five spins Si^ ... Sm,, are on the five vertices of pyra¬ 
mid /r. The parameter J sets the energy scale, and the 
linear system size is L, with a total of N = spins. 
Pyramids pi for which the product of spins is equal to 
— 1 are defeets which carry an energy J: from a ther¬ 
modynamic viewpoint they are non-interacting and their 
statistics are those of an ideal gas. Hence for a single sys¬ 
tem there are no thermodynamic phase transitions at any 
finite temperature. However, the dynamical behavior of 
the defects is complex and co-operative: by analogy with 
a similar model in two dimensions m, it is believed that 
for low temperatures T, the relaxation time r diverges 
Inr ^ T~^ [16] . and the range of certain multi-point 
spin correlations also diverges [31] . 

The overlap is Q{C,C') = which measures 

the similarity between spin configurations. We draw a 
configuration C' at random from an equilibrium distri¬ 
bution at temperature T'. Holding C fixed, we calcu¬ 
late expectation values with respect to the distribution 
p{C\C') = Z{C')-^ Q[eNQ(Cfi')-Eo{C)]/T ^ ^{C') is a 

partition function. Thus for e = 0 we consider equilib¬ 
rium behaviour of the isolated SPyM at temperature T, 
while increasing the coupling e biases the configuration C 
to increase its overlap with the reference C. Finally we 
perform an average over the reference configuration C. 
We set J = 1 which fixes our energy unit, so the dimen¬ 
sionless parameters of the system are (T, T',e). We use 
a Monte Carlo sampling scheme to study these coupled 
systems, as described in Appendix A. 

Results -Fig.j^shows results for T' = T, in which case 
the reference configuration is representative of thermal 
equilibrium. The mean overlap {Q) in Fig. j^a) increases 
with e, as expected. This increase is gradual for small e, 
before a steep increase sets in at larger e. The theoretical 


2 


prediction [illllllH] is that the gradient d{Q)/de should 
diverge at an RFIM critical point at some {Tc,ec), and 
that for T < T^, one should observe (in the thermody¬ 
namic limit) a first-order transition, with a discontinuous 
jump in (Q). In these finite systems, no divergences are 
observed: to demonstrate the existence of the phase tran¬ 
sition we use finite-size scaling methods. 

The reference configuration C is a source of quenched 
disorder in this problem, and averages are calculated in 
two stages: first a thermal average at fixed C, which we 
denote by {■)c> and then a disorder average, denoted by 
(•). The notation (•) = (•)c' indicates the double average. 
Fig[^b) shows the average (Q) as well as the behaviour 
of (Q)c' for eight representative configurations C'. For 
each C', one sees a very sharp jump in the overlap at a 
sample-dependent field However, the jump in the 
average overlap is broadened out due to fluctuations in 
ej/ between samples. This is the expected behaviour at 
an RFIM critical point [22] “ it means that the disorder 
is the dominant source of fluctuations in the problem. 

Fig. [^c) shows the (total) susceptibility Xtot = 
N{SQ'^) = N[{Q'^) — (Q)]. As the temperature is reduced, 
the fluctuations grow rapidly and depend increasingly 
strongly on the system size, as expected in the vicinity of 
a phase transition. Finally, Fig. m shows the distribu¬ 
tion of the overlap Q, evaluated at the field s* which max¬ 
imises Xtot- For the lower temperatures, a two-peaked 
structure is clearly visible. Considering these data to¬ 
gether with Fig. Sb), one sees that typical reference 
configurations contribute to either the low-Q or high-Q 
peak, as expected for RFIM criticality: the proportion of 
configurations C' that contribute simultaneously to both 
peaks scales as L~^ [22| where d is a critical exponent for 
the RFIM (in three dimensions, 9 « 1.5 [331|21|). Also, 
P{Q) in Fig. j^d) has an asymmetric shape, with a very 
broad high-Q peak. From Fig. , we attribute this 
broadness to the fluctuations between reference config¬ 
urations: the position and the size of the jump in {Q) 
vary signficantly. Finally, we note that numerical uncer¬ 
tainties in P{Q) are considerable: where the probability 
is small our estimate of P{Q) may be dominated by just 
one or two out of the 64 reference configurations consid¬ 
ered. We return to this point below. 

Taken together, the data in Fig. are consistent with 
an RFIM critical point for T' = T in the range 0.3 — 0.4. 
However, while the MC methods used here lead to ef¬ 
ficient simulations even close to phase transitions (see 
Appendix A), the slow (glassy) relaxation of the model 
limits the temperatures and system sizes that we can 
consider. If an RFIM critical point exists for T' = T, we 
expect similar critical points at other reference temper¬ 
atures T'. We therefore analyse a lower reference tem¬ 
perature T' = 0.25, providing further evidence for RFIM 
criticality in this system. 

For T' = 0.25, Fig.j^a) shows a sharp jump in {Q) as e 
increases, this time at a relatively high temperature, T « 
0.8, for which simulations are more tractable. This allows 
investigations of larger systems, up to fV = 32^. Fig.j^b) 



FIG. 1. Results for T = T'. (a) Mean overlap (Q) as a func¬ 
tion of £, showing an increasingly sharp jump as T is reduced. 
The dashed line is {Q) = e/T, which gives the linear response 
behaviour around e = 0. The system size is = 16®. Colors 
and symbols have the same meaning in all panels, (b) Average 
overlap for T = 0.333, as well as thermal averages {Q)c' for 
eight representative reference configurations, (c) Total sus¬ 
ceptibility Xtot = N{5Q‘^). Solid symbols are for a system of 
size N = 16® while open symbols are for = 8®, at the same 
temperatures, (d) Distributions P{Q) for N = 16®, at the val¬ 
ues of e which maximise Xtot, showing increased bimodality 
on reducing temperature. 

shows the distribution of the overlap P{Q), whose form 
changes from unimodal to bimodal as the temperature is 
lowered. 

To identify the critical point, we use finite-size scal¬ 
ing. However, the RFIM has unusual finite-size scaling 
properties: the transition is almost first order and the 
distribution P{Q) has two well-separated peaks even at 
the critical point [23] . The order parameter exponent /3 
is also very close to zero [24]. Together with the signifi¬ 
cant numerical uncertainties in P{Q), these two features 
lead to difficulties with classical finite-size scaling based 
on universal cumulant ratios or order parameter distribu¬ 
tions. To address this problem, we follow [23]: for each 
reference configuration C', we calculate the coupling 
that maximises the variance {SQ^)c'- At this coupling, 
the distribution of Q is typically bimodal, even above 
Tc [23]. We define a free energy for this reference con¬ 
figuration P*{Q) = —TlnPc'{Q). This free energy has 
two minima, and the height of the barrier between them 
is the interfacial free energy cost between high- and low- 
overlap states. To obtain an average free energy cost, we 
calculate the average free energy P*{Q) from which we 
extract the barrier height AP. 

Fig-d c) shows the resulting free energy (scaled by 
/? = 1/T to allow interpretation as a log-probability). In 
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FIG. 3. (a) Free energy j3J-*{Q) for T = T' = 0.333 and 
system sizes N = (8®, 16^). The interfacial cost increases 
strongly with system size. The system is at coexistence; the 
difference in height between the two minima reflects their 
different widths, (b) Scaled interfacial costs filled 

symbols are N — 16^ and open symbols N = 8^. The costs 
are non-monotonic in T; their finite-size scaling is discussed 
in the main text. 


FIG. 2. Results for T' = 0.25. (a) The mean overlap shows a 
sharp crossover as e is increased. The system size is Ai = 32®. 
The dashed line is the linear response relation {Q) = e/T. 

(b) Overlap distribution P{Q), evaluated at e = e*; the sym¬ 
bols and colors indicate the same temperatures shown in (a). 

(c) Free energy I3J-*{Q) for T — 0.833, showing that the in¬ 
terfacial free energy cost (barrier between the two minima) 
increases strongly with system size. Error bars show standard 
errors for N = 32®; numerical uncertainties for N = 16® are 
comparable with system sizes, (d) Temperature dependence 
of the scaled interfacial free energy cost, /3AF’/L®, from which 
we estimate Tc « 0.83. 


contrast to P{Q), the numerical uncertainties in J-*(Q) 
are straightforward to estimate, and relatively small. For 
RFIM critical points, one expects a sigificant interfa¬ 
cial free energy in small systems even for T > Tc, but 
lim 2 ,_>oo 0 in that case. For T < Tc, one expects a 

divergence of the interfacial cost as ~ as usual 

for first order phase transitions (here, d = 3 is the spatial 
dimension). At the RFIM critical point, the interfacial 
free energy cost AT scales as A® [^[13]. In Fig. j^d), 
we plot PATjL^ for two system sizes (we take 9 = 1.5). 
The interfacial costs are large and they grow with sys¬ 
tem size, Fig. [^b); they grow increasingly rapidly at low 
temperatures, Fig. [^c). This presents strong evidence 
for the existence of a critical point, whose properties are 
consistent with RFIM universality. From the crossing 
point in Fig. [^d), we estimate the critical temperature 
as Tc ~ 0.83. 

Interfacial costs for T' = T. Armed with these re¬ 
sults for T' = 0.25, we now return to the important 
case T' = T. Fig. [^a) shows the free energy T*{Q), at 
T = T' = 0.333. The interfacial cost increases strongly 
with system size. Fig. Sb) shows the temperature- 
dependence of this cost. For T' = 0.357 and T > T', 
the scaled free energy cost always decreases as the sys¬ 
tem size is increased from N = 8^ to N = 16®. This 


indicates that these systems are all above Tc- We also 
note that these free energy costs are non-monotonic in 
temperature - our interpretation of this unusual feature 
is that the high-overlap and low-overlap states have in¬ 
creasingly similar structures as T approaches T', which 
tends to reduce the surface tension between them. More¬ 
over, for T' = 0.333, the scaled free energy costs for the 
two system sizes are very close to each other over a wide 
range of temperature T. [This is in contrast to the clear 
crossing of the curves in Fig. [^d).] Our interpretation 
of this last result is that all these systems are close to 
criticality. 

Summary of phase behaviour - To illustrate the result¬ 
ing scenario we consider the critical temperature Tc as a 
function of the reference temperature T' . For T' = 0 the 
reference configuration C is in its ground state which has 
all spins up, and the system reduces to a SPyM in a mag¬ 
netic field, for which there is known to be an Ising critical 
point at Tc = 0.98 [Ull^. On increasing T’, we expect a 
line of critical points in the parameter space (e, T, T'), all 
of which should be of RFIM type, except for the special 
(Ising) case T' = 0. The estimated critical ponts that we 
have found in this study are shown in Fig. |4];a). They 
separate a region of parameter space in which phase coex¬ 
istence is possible (“two phase”) from a one-phase region 
where the response to the coupling e is smooth through 
the entire range of the overlap (0 < Q < 1). The Franz- 
Parisi potential V(Q) is strictly convex in the one-phase 
region but includes a linear segment (Maxwell construc¬ 
tion) in the two-phase region. (Evidence for the critical 
point with T' = 0.286 is shown in Appendix B.) 

The dashed line in Fig. |^a) indicates T = T'. It inter¬ 
sects the line of critical points at T' « 0.33. When coex¬ 
istence takes place at T = T', the nature of the interfacial 
free energy cost is a subtle one, which has important im¬ 
plications for theories of the glass transition. Since C is 
an equilibrium configuration of the isolated system, one 
expects the high-Q state to have configurations C that are 
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FIG. 4. (a) Estimated critical temperatures as a func¬ 

tion of T'. These critical points delimit a two-phase region 
(shaded) within which the Franz-Parisi potential V{Q) in¬ 
cludes a Maxwell construction, associated with phase coex¬ 
istence. The line of critical points intersects the dashed line 
T — T' at T ~ 0.33. (b) Sketched phase diagram for the 
SPyM with T' = T. The points indicate the susceptibil¬ 
ity maxima from Fig. with the estimated critical point at 
T « 0.33. Lines follow e = (see main text) with fit 

parameter A = 1.25. They indicate a first-order phase tran¬ 
sition (dashed) and a susceptibility maximum (dotted). 


also close to equilibrium. Similarly, the low-Q states are 
only weakly perturbed from equilibrium, so one expects 
C to have close-to-equilibrium structures in that case too. 
Thus, the phase coexistence that is illustrated in Fig.j^is 
taking place between different states which have (statis¬ 
tically) very similar structures. The interfacial cost as¬ 
sociated with such states plays a central role in random 
first order transition (RFOT) theory [26], because C may 
be viewed as a form of self-generated disorder. However, 
the nature of these interfaces in finite-dimensional sys¬ 
tems is not fully understood |27||28]. As noted above, 
the reduction in interfacial cost as T approaches T' in 
Fig. I^b) is consistent with the idea that the high- and 
low-overlap states are becoming increasingly similar. It 
also indicates that the cost is primarily entropic and not 
energetic in origin (otherwise the probability cost would 
likely increase on cooling, as 13/S.F = AE/T). 

Still concentrating on the case T = T', the phase co¬ 
existence for T < Tc reveals that the system supports 
many metastable states, with a finite interfacial tension 
between them. However, we emphasize that for e = 0 
the free energy of the SPyM is a simple analytic func¬ 
tion of temperature and there are no finite-temperature 
phase transitions. In general one expects a first-order 
transition at e*/T « ScAQ where Sc is the configura¬ 
tional entropy density and AQ the jump in the overlap 
at £* [29]. For plaquette models one expects Sc to be of 
the same order as the total entropy s = e“^/^(l/r — 1), 
and AQ 1, which gives an estimate of the critical 

coupling, e* ~ In Fig. |^b), we sketch the ex¬ 

pected phase behavior of the SPyM, including three data 
points that indicate the locus of e* as T is varied, hold¬ 


ing T’ = T (recall Fig. ^a))- Given the limited data, 
the fitted line e* ^ e“^/^ should be regarded only as a 
qualitative prediction for the low temperature behaviour, 
but we emphasise that the first-order line can meet the 
£ = 0 axis only at T = 0. (Similar scaling is found for an¬ 
nealed coupling between plaquette models, although the 
phase transitions are of Ising type (not RFIM) in that 
case [TFlIBP] .) 

Discussion - The behaviour of the SPyM is character¬ 
istic of the DF perspective of the glass transition (ab¬ 
sence of a thermodynamic transition, effective kinetic 
constraints, facilitated relaxation and heterogeneous dy¬ 
namics). We have shown, nevertheless, that when cou¬ 
pled to a quenched reference state, the SPyM displays 
thermodynamic overlap transitions in the RFIM class, 
as expected from mean-field calculations and RFOT the¬ 
ory. It follows that quenched overlap transitions in atom¬ 
istic simulations m are not incompatible with DF the¬ 
ory; also transitions occurring at T, £ > 0 do not im¬ 
ply any finite-temperature singularities at £ = 0 (recall 
Fig.|^b)). More generally, and taken together with other 
recent studies of plaquette models [5TH55] . the presence 
of overlap transitions in the SPyM shows how these mod¬ 
els can be consistent at the same time with predictions of 
both DF and RFOT theory, offering a link between these 
scenarios. 

To forge such a link, the crucial question is whether 
(and how) phase transitions for £, T > 0 are related to 
the unbiased {e = 0) properties of the system. Given 
the first-order transition line in the e-T plane, we ex¬ 
pect a nucleation-like mechanism for dynamical relax¬ 
ation near this line, with an initially high-overlap state 
decaying to a low-overlap one via nucleation and growth 
of a droplet of the low-overlap phase. This situation 
is close to the RFOT picture described in [50]: for 
£ < £* a nucleation argument yields a free energy barrier 
PAF ~ (£* — £)i?c ~ (£* — ey~‘^ where Rc ^ l/(£* — c) 
is the size of the critical nucleus, which diverges at the 
first order transition line. However, for £ = 0 the nat¬ 
ural picture is that associated with DF: point-like de¬ 
fects facilitate co-operative rearrangements over a length 
scale ^ ~ where d{ is a scaling exponent |14j . 

The free energy barrier for these processes scales as 
PAF ~ (l/r)ln^ ^ (1/T)^. Bridging between the two 
regimes £ ~ £* and £ ~ 0 remains a challenge. However, 
our results suggest that such an interpolation might be 
a way to connect the defect-mediated dynamics of DF 
theory to the nucleation picture predicted by RFOT and 
mean-held theories. 
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Appendix A: Numerical methods 


1. Thermal averages 


Our methods for simulating the SPyM follow closely 
those of [16]. We start by using the continuous time 
Monte Carlo (CTMC) method of [M] to simulate the 
model, using Glauber rates for spin flips that take into 
account only the uncoupled energy Eq. We modify this 
algorithm slightly, in order to sample ensembles where 
the distribution of the configuration C is 


Pb(C|C') OC e-^o(e)/T+eNQiC,C')/T 


1 

b{Q{C,C'))- 


(Al) 


Here b{Q) is a bias function that is chosen to promote 
rapid sampling of a wide range of Q-values in a single 
simulation [35j . The construction of a suitable choice 

for b is discussed below. _ 

To sample the distribution ( |Al[ ) within the CTMC dy¬ 
namics, we incorporate an additional accept/reject step 
for each move, using a Metropolis criterion 


Tacc = min < 1, exp [{eN/T){Q 

new ^old)] 


^(Qold) 

b{Q 

new) 

(A2) 

It is easily verified that this Monte Carlo algorithm con¬ 
verges to a steady state whose a distribution is given by 
(Al). The desired Boltzmann distribution is then ob¬ 
tained by a reweighting p(fi\C') oc ph{C\C') ■ b{Q{C,C')). 

We use an automated algorithm for choosing and up¬ 
dating the bias b{Q). The method operates for a fixed 
reference configuration C', so the final (quenched) aver¬ 
age over C' simply involves running the algorithm many 
times for a range of different reference configurations. 

We focus on the biased distribution of the overlap 
PhiQ\C') that is obtained from pb(C|C'). Clearly, the 
biased distribution is Ph{Q\C') oc Pc'{Q)/b{Q) where 
Pc'iQ) is the unbiased overlap distribution discussed in 
the main text. For systems close to first-order phase tran¬ 
sitions, there is a long timescale associated with sponta¬ 
neous transitions between the phases. The time required 
for such a transition to occur spontaneously scales with 
the depth of the minimum in Pc {Q) and makes it difficult 
to assess the precise location of the phase transition (due 
to hysteresis effects). The purpose of the bias function 
b{Q) is to smooth out any deep minima in P{Q), so that 
the distribution Ph{Q) is almost flat over the range of Q 
that is of interest [35]. That is we choose, b{Q) « Pc'iQ)- 
In practice, we start at a relatively high temperature 
T for which sampling is easy and we estimate the value 
of e such that Pc{Q) has maximal variance. We use the 
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MC dynamics with 6 = 1 to collect Ns representative 
configurations C'' with r = 1,2,..., Ng. We take Ng in 
the range 10^ — 10"^: samples are not all independent 
from each other but the sampling runs are long enough 
that the system fully decorrelates within each run. For 
each sample, we store the overlap Qr = Q(C^,C') and 
the energy = Eo(C’'). This provides an estimate for 
a suitable bias potential for further simulations at this 
state point: 


b(Q) (xJ^Sq^q^ (A3) 

r 


where (5 q,q' is the Kronecker delta. Hence b{Q) is the 
empirical distribution of Q. In practice, Eq ( |A3[ ) can be 
used only in the range of Q for which reliable data are 
available, and we take b = const, outside this range. 

Now reduce the temperature to T — AT. We also re¬ 
duce the coupling e to a value e' whose choice will be 
described below. We can estimate Pc'{Q) at this new 
state point as 


Pc'{Q) « b\{Q) oc 


exp(-AjA^ AgNQr) 


(A4) 

where Aj = 1/(T — AT) — 1/T is the change in inverse 
temperature and A^ = s' jlT — AT) — e/T. The param¬ 
eter e' (or equivalently A^) is chosen to maximise the 
variance of b/\(Q), so that the new state point remains 
close to the susceptibility maximum as the temperature 
is reduced. The parameter AT must be chosen small 
enough that &a(Q) provides a reasonably accurate esti¬ 
mate of Pc'iQ) at the new state point: in practice this 
means that the exponential weights in (A31 should not 
lead to to a sum dominated by just a few terms. 

Repeating this procedure allows the temperature to be 
reduced while b{Q) is maintained close to Pc'{Q)- This 
is the method used to obtain the strongly bimodal dis¬ 
tributions P{Q) in Figs Id, 2b, etc. of the main text. 

We note in passing that the Franz-Parisi potential may 
be obtained as 


V{Q)=T \im N-hnPc'{Q). (A5) 

Af—>00 

This function is therefore directly accessible by our 
method. It is related to the free energy T*{Q) of the 
main text as 


V{Q)= lim 

N^oo 




s*Q 


2. Averages over reference configurations 

Plaquette models such as the SPyM have the useful 
feature that equilibrated reference configurations can be 
generated in a simple way by choosing random defect po¬ 
sitions at a density l/(l-|-e^'^'^) and then constructing the 
relevant spin configuration as described in [16]. The pro¬ 
cedure described in the previous section is then repeated 
for each reference configuration, leading to the averaged 
quantities shown in the main text. We typically calcu¬ 
lated averages over 64 reference configurations, except in 
cases where fewer samples were sufficient to obtain con¬ 
vergence. 


3. Numerical parameters 

In the main text we quote all temperature to three 
significant figures. The exact temperatures used were 
T = 1/2.6,1/2.8,1/3.0 (Fig. 1); T = 1/1.1,1/1.2 and 
T' = 1/4 (Fig. 2). In Fig. 4 there is a single data point 
with T = 1/3.5 « 0.286. 


Appendix B: Results for T = 0.286 

Fig. [5] shows the scaled interfacial cost in a system 
with reference temperature T' = 0.286, comparing sys¬ 
tem sizes N = 16^,32^. Compare Fig. 2d of the main 
text. We infer that the critical temperature in this case 
is Tc « 0.71: the position of this critical point is shown 
in Fig. 4a of the main text. 


T’ = 0.286 



FIG. 5. Results for the scaled interfacial cost as a function 
of temperature T, with reference temperature T' = 0.286. 
Compare Fig. 2d of the main text. 
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(A6) 









